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Abstract. Here we introduce the concept of "optimal particles" for strong 
interactions with electromagnetic fields. We assume that a particle occupies a given 
electrically small volume in space and study the required optimal relations between the 
particle polarizabilities. In these optimal particles, the inclusion shape and material are 
chosen so that the particles extract the maximum possible power from given incident 
fields. It appears that for different excitation scenarios the optimal particles are 
bianisotropic chiral, omega, moving, and Tellegen particles. The optimal dimensions 
of resonant canonical chiral and omega particles are found analytically. Such optimal 
particles have extreme properties in scattering (for example, zero backscattering or 
invisibility) . Planar arrays of optimal particles possess extreme properties in reflection 
and transmission (e.g., total absorption or magnetic- wall response), and volumetric 
composites of optimal particles realize, for example, such extreme materials as the 
chiral nihility medium. 
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1. Introduction 

Let us consider an electrically small particle excited by a given external electromagnetic 
field, for example, a plane wave. We assume that the particle is electrically small and 
its response can be modelled in the dipole (electric and magnetic) approximation. The 
incident fields polarize the particle, and we are interested in the question how this effect 
can be maximized. Of course, simply increasing the values of the particle polarizabilities 
(for example, increasing the particle size or bringing it to a resonance), we increase the 
amplitudes of the excited dipole moments. But what would be the optimal relation 
between the polarizabilities of the particle to ensure the most effective interaction with 
the incident fields, provided that the overall size of the particle is fixed and we cannot 
increase the absolute values of the polarizabilities? 

Recently, it was noticed that the shape of chiral particles (for example, metal 
helices) can be optimized in such a way that these optimal helices radiate waves of 
only one circular polarization, whatever is the polarization of the exciting field [H [2] . 
Later, it was shown that the optimal spiral shape corresponds to the maximum (or 
minimum, equal zero) reactive energy of the particle in a given plane- wave field [3], 
interacting with the wave in the most effective way. It has been therefore established that 
there exist particles (and metamaterials formed by these particles), whose properties 
are optimal for strong interactions electromagnetic fields. Furthermore, mixtures of 
optimal spirals realize effective media with the index of refraction n = — 1 for one of the 
circular polarizations, while the same medium is transparent for the orthogonal circular 
polarization [1]. Analytical and numerical studies of optimal helical particles have 
been done in [5]. However, optimal helices show their optimal and extreme properties 
only if excited by propagating circularly polarized plane waves. First steps towards 
generalizations for other exciting fields have been done in [6l [7], where the case of 
linearly polarized evanescent fields has been discussed. 

The goal of the present work is to find the optimal particles which would interact 
most effectively with waves of other polarizations (in particular, linear) and also with 
both propagating and evanescent incident fields. We expect that such optimal inclusions 
will open a way to realizing particle arrays and volumetric composites with extreme 
properties, as we see from the known example of the optimal spiral. 

The strength of field-particle interactions can be measured in various ways. Here 
we define the optimal particle as one which extracts the maximum (or minimum) power 
from given incident fields. 

2. General framework 

From the known example of the optimal chiral particle it is clear that in order to enhance 
the field-particle interaction, we need to design a particle which would be polarizable 
by both electric and magnetic fields, and the particle should be bianisotropic, and not 
just an electrically and magnetically polarizable inclusion. In the presence of magneto- 
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electric coupling, excitations of the particle due to the incident electric and magnetic 
fields can sum up coherently, maximizing the total excitation. Thus, we assume that the 
particle for optimally strong interactions with different kinds of waves is a bianisotropic 
particle whose electric and magnetic moments are induced by both electric and magnetic 
fields E and H. The general polarizability relation can be written as 
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where p and m are the induced electric and magnetic dipole moments, respectively, and 
Oij are the polarizability dyadics. 

Let us define a direction is space (along the unit vector zq) and assume that the 
particle is excited by a plane wave propagating along this direction. In other scenarios, 
we can assume that the particles form a planar array in the plane orthogonal to zq, or 
that the particles form a 3D periodical crystal forming a half space (or a finite-thickness 
layer), with the surface orthogonal to zq. In all these scenarios, the most interesting 
case is the case of the uniaxial symmetry of the particles and the polarizability dyadics. 
To simplify the following derivations, we will in addition assume that the particle is not 
excited by the fields directed along the axis zq. This fully covers the case of excitation 
by plane waves propagating along the axis zq and the case of the normal incidence on 
planar arrays and planar interfaces. For the general excitation, the assumption means 
that the theory is restricted to "planar" particles, which feel only the fields in the plane, 
orthogonal to zq. Thus, we write the polarizability dyadics of uniaxial bianisotropic 
particles in the form 

^— CO ~r I cr T V / 

Here, It = /— zqZq is the transverse unit dyadic, I is the 3D unit dyadics, and = zq x 
is the vector-product operator. 

In the following, classification of bianisotropic particles in terms of reciprocity and 
the symmetry of their magneto-electric coefficients will be important. We will discuss 
the two reciprocal classes (chiral and omega) and two nonreciprocal classes (moving 
and Tellegen) [8]. For reciprocal particles the electric and magnetic polarizabilities are 
always symmetric dyadics, so we keep only a™ and a'^^ in the expressions for clee 
and amm- For the nonreciprocal classes we keep also the anti-symmetric parts aH and 
•^mm of these dyadics. The polarizability dyadics of particles of these four classes are 
summarized in Table [H 

The imaginary units in the notations for the coupling parameters of reciprocal particles 
and in the anti-symmetric parts of electric and magnetic polarizabilities are introduced in 
order to ensure that the absorption and scattering losses are measured by the imaginary 
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Table 1. Uniaxial bianisotropic particles 



Main classes of uniaxial bianisotropic particles 


Omega 


Chiral 


Moving 


Tellegen 
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parts of all the polarizability components. This follows from the general conditions for 
lossless particles (e.g. [8]): 

Q^mm '-'^mm' (^) 

where f is the Hermitian conjugate operator (transpose of the complex conjugate). 

Next, we specify the assumptions about incident electromagnetic fields. We 
will study excitations in form of plane waves (propagating or evanescent) or their 
combinations (standing waves). As the particle reacts on the incident electric and 
magnetic fields in the x — y plane, we separate the fields of a general plane wave in free 
space into transverse and normal field components: 

H = Ht + -ff^zo, 

(4) 

E = Et + EzZq, 

where 

Ht = -ff^xo + HyYo, 

(5) 

Et = -Z • zo X Ht. 

The wave vector k of a plane wave we also split into the longitudinal part /3zo and 
the transverse part kt: k = kt + /3zo, where /3 = (cj^/iQeo — k'^Y^'^ . The dependence 
on the coordinates is exp(— jk • r), and the harmonic time dependence is in the form 
exp(ja;t). The wave impedance relates the transverse components of the fields, and it 
can be written as 



ktkt , ^ kt X zokt X zq 

/3 
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(6) 
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Here, Ztm and Zte are the characteristic impedances for the two linear 
eigenpolarizations (transverse electric, TE, and transverse magnetic, TM, waves are 
defined with respect to the vector zq). 
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Next, we write the general expression for the time- averaged power extracted by the 
particle from the incident field (the external field power spent on the excitation of the 
particle): 

P = -|lm{p-E: + m-H:}, (7) 

where * denotes the complex conjugate operator. Using ([1]) and ([2]), the electric and 
magnetic moments of bianisotropic particles can be generally written as 

m = a^^^Et + a^^zo x Et + a^^^Ht + a^mZo x Hf 

We express the transverse components of the incident field in the Cartesian coordinates 
X, y, and substitute these moments in ([7]). The result reads 



P = --lm{a::\E, 



+a^' ( ^—\E P ^—\E P 

+r.co I J—R R* J —E*E 1 (9) 
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This is the power absorbed by a bianisotropic particle in a given arbitrary polarized field. 
In what follows, we will find the optimal relations between the particle polarizabilities 
which maximize this value, providing recipes for the design of particles for extremely 
strong interactions with electromagnetic fields. 



3. Optimal particle classes 

Let us assume that the amplitudes of the polarizabilities in ([9]) are limited or fixed by the 
particle size and material. In this situation, it is clear that the topology of the optimal 
particle should be chosen so that the terms containing the magneto-electric coefficients 
would contribute constructively to the total extracted power. This consideration will 
determine the optimal particle classes for different types of excitation, because the 
impedance values in ([9]) can take real or imaginary values for propagating waves or 
reactive external fields. 
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3.1. Excitation by linearly polarized propagating waves 

If the particle is excited by a linearly polarized propagating plane wave, then the two 
transverse field components and Ey are in phase or out of phase {E^E* is a real 
number), and the wave impedances Zte and Ztm are real numbers. The expression for 
the absorbed power in ([9]) takes the form 

|2t r„,co 



P = --|Et|^ImK: 
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The terms related to omega and chiral coupling do not contribute to ffTOl) and the 
remaining co- and cross-coupling polarizabilities refer to the Tellegen and moving 
particles, respectively. It means that for absorbing the maximally possible power from 
incident propagating linearly polarized waves, particles with nonreciprocal magneto- 
electric coupling are required (both Tellegen and "moving" coupling mechanisms). In 
the case of an axially incident wave (propagating along zq) both impedances are equal to 
r/o = the term concerning to the Tellegen coupling cancels out, and the absorbed 
power can be written as 

P = -^|E,pIm |< --V + ^a-^l . (11) 

Therefore, in the case of an axially incident linearly polarized propagating wave, the 
optimal particle class is that of moving particles (nonreciprocal anti-symmetric magneto- 
electric coupling). 

3.2. Excitation by linearly polarized reactive fields 

In case of reactive linearly polarized field in free space (for example, an evanescent plane 
wave or a standing linearly polarized wave) the transverse electric and magnetic fields 
are 90° out of phase and both impedances Zte and Ztm are purely imaginary. The 
absorbed power in ([9]) can be written as 

'2t r„,co 
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We see that in this case the absorbed power can be enhanced introducing reciprocal 
chiral and omega couphng effects. If the two impedances are equal Zte = -^tm = Z 
(for example, a standing plane wave at the plane equally distanced from the maxima of 
the two fields), the term with the chirality parameter cancels out and we can write an 
expression analogous to that for the propagating waves (ITTl) : 

P = -||E,pIm |< + 2mm (i) + • (13) 

We can observe that for optimizing interactions with propagating linearly polarized 
waves, nonreciprocal particles are needed, whereas in the case of linearly polarized 
reactive fields, we need reciprocal (omega) particles. It is interesting to note that the 
omega particle was introduced as a particle strongly interacting with linearly polarized 
fields P [To] , but it was noticed that they are optimal for linearly polarized evanescent 
fields. 



3. 3. Excitation by circularly polarized propagating waves 

In the general case, we define the circularly polarized incident field as a field 
configuration where the electric field vector is circularly polarized in the transverse 
plane of the particle, that is, Ey = ^jE^ [E^E* is an imaginary number). For circularly 
polarized waves, the upper and lower signs correspond to left- and right-hand circularly 
polarized waves propagating along zq, respectively. In general, magnetic field can be 
elliptically polarized. 

For propagating waves, the impedances .^te.tm are real numbers, and the power 
extracted by the particle according to ([9]) can be written as 



P = --|Etrim{a 



ee 




(14) 



Here, the contributing co- and cross-coupling polarizabilities stand for reciprocal (chiral) 
and nonreciprocal (moving) classes, and the terms concerning to omega and Tellegen 
properties do not have any role in interacting with circularly polarized propagating 
waves. In addition, the antisymmetric parts of electric and magnetic polarizabilities 
contribute to the received power. As it can be seen from ( IT^ . in the case of axially 
incident waves, both chiral and moving-particle coupling terms are essential, and the 
power can be written as 

P = -||E,pIm |< - ^ (y T «:) + ± U: + ^) } • (15) 
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This is in agreement with the earher studies of chiral particles as the optimal particles for 
interactions with propagating circularly polarized waves [H El H] , where it was assumed 
that the particles were reciprocal (that is, parameters V, aH and a"^^ equal zero). 

3.4- Excitation by circularly polarized evanescent fields 

In this case, we assume that the wave impedances are purely imaginary and find that 
the absorbed energy in ([9]) can be simplified to 

^ |2t™ r CO 



P = --|Et|'lm{« 
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Similarly to the case of propagating circularly polarized waves, both reciprocal and 
nonreciprocal magneto-electric coupling contribute to the extracted power, but here we 
need the other two classes of magneto-electric coupling: reciprocal omega coupling and 
nonreciprocal Tellegen coupling. 



4. Optimal reciprocal particles 

In the previous section, we found what classes of magneto-electric particles are most 
suitable for interactions with different types of exciting fields (propagating versus 
evanescent, linearly versus circularly polarized). The next step is to find the optimal 
particle shapes for these particles. In this section, we will do it for reciprocal particles: 
chiral and omega particles, optimal for circularly polarized propagating waves and 
evanescent linearly polarized fields. Optimal chiral particles have been studied earlier 
[U m El m |5], but in those papers the adopted criteria for the optimal field-particle 
interaction were different: Earlier introduced optimal helices react only on circularly 
polarized waves of one of the two orthogonal polarizations and store the maximum 
reactive field defined in the quasi-static approximation. With this in view, we will 
compare our results for the optimal dimensions of helices optimized in the sense of the 
maximum extracted power from the incident fields. 

The optimal shape of electrically small reciprocal particles is possible to determine 
under quite general assumptions thanks to the existence of so called "hierarchy 
of polarizabilities" . Considering the particle size (d) normalized to the free-space 
wavelength d/X as a small parameter, it is known that the electric polarizability 
is the strongest (exists even in the zero-frequency limit), then comes the cross- 
coupling coefficient (a first-order spatial dispersion effect of the order of d/\), and the 
weakest effect is the artificial magnetism (a second-order effect of the order of (d/A)^) 
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[TT| [121 [ini El El- To estimate the polarizabilities and find the optimal particle shape 
analytically we will use the antenna model for thin-wire scatterers (chiral and omega) 
[H]. We will see that the solution for the optimal particle shape will define an optimal 
value for the "scaling parameter" for the polarizabilities of reciprocal particles. 



4.I. Chiral particles 

According to the above results, chiral particles can be optimized for interactions with 
propagating circularly polarized waves. As a particular example of a small chiral particle, 
interacting with fields in the microwave region, we will consider the well-known metal 
canonical particle shown in Figure [1] This particle can be viewed as a connection of 
two small wire antennas: a short electric dipole antenna of length 2/ ^ A and a small 
loop antenna with the radius a ^ A [14J. The plane of the loop is orthogonal to the 
electric dipole antenna. This simple shape of a metal helix was probably introduced 
for the first time by D. Jaggard et al. [15] and then used in a number of studies as a 
canonical example of a chiral particle which allows simple analytical solutions for its 
polarizabilities (see e.g. [HI E] and references therein). 




Figure 1. Canonical chiral particle. A uniaxial particle with the polarizabilities 
defined in Table [T] can be formed by adding a second identical particle rotated by 90° 
around the axis z. 



In terms of the parameters of the two small antennas forming the particle, the 
polarizabilities along axis x can be written as 
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Here, Z^ot is the total impedance of the particle (the sum of the input impedances of 
the loop and electric dipole antennas), 5* = vra^ is the loop area, and = u^/e^JIo is 
the free-space wave number. At the resonance of the particle, Z^^^t = -Rtot + j-^tot is a 
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real number, which is non-zero because of inevitably present scattering losses even in 
the absence of absorption in the particle. Thus, at the resonance frequency, the electric 
and magnetic polarizabilities Oee, «mm, and the coupling (chirality) parameter k are 
imaginary and, as it can be seen from (fT^ . all of them contribute to absorption of field 
energy. 

The absorbed power f|T^ can be written as 



Substituting the characteristic impedances from ([6]), we get the absorbed power at the 
resonance frequency as 

We simplify our problem assuming that the incidence direction is orthogonal to the 
induced electric and magnetic moments (along the axis z of Figure [1]), which means 
(3 = ko- The absorbed power simplifies to 

2 Rtot I I \ I 



P=Ue,\'^{i±2^+^\ . (20) 



For particles without absorption loss, i?tot is only due to the scattering loss, and its 
value is known from the theory of small antennas (e.g., [TOl ITT]): 

Substituting this value, after some simple algebra, we find P as 

2 



Rtot = ^ {kll' + ktS') . (21) 
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Here we have used the resonance condition for the particle: at the fundamental resonance 
frequency the total (stretched) length 21 + 27ra of the wire is approximately A/2. To 
determine the optimal relation between / and a which would bring the received power 
P to an extremum, we calculate its derivative with respect to / and equate it to zero 
under the condition of resonance. This determines the optimal dimensions: 

lopt = (2 - V3)^, 

(23) 

aopt- ^ ^. 

This corresponds to the optimal value of the "scaling factor" k^S/l = ±1 in (120!) . which 

means that all three normalized polarizabilities of optimal chiral particles are equal: 

1 , , 

rjoaee = =F/t = — amm- (24) 
Vo 
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This last result coincides with that in [H [31 S] and the optimal dimensions f l23p are 
the same as found in [5J, but an all these earlier papers the optimal particles were 
defined as those interacting only with circularly polarized waves of one of the two 
eigenpolarizations. Now, we see that these particles are also optimal in the sense of 
the maximum received power. The physical nature of this optimal behaviour of helices 
is simple: In one extreme case the electromotive forces induced by electric and magnetic 
incident fields cancel out, and the total induced current is zero, while in the other 
extreme case the two forces sum up in phase and interaction is optimally enhanced. 
Thus, we see that all the known definitions of optimal chiral particles lead to the same 
solution for the particle shape and to the same general balance relation between the 
polarizabilities (El]), if the incident wave propagates along the particle axis. Note that 
in general, the optimal particle dimensions depend not only on the frequency, but also 
on the transverse wave number (the incidence direction) of the exciting plane wave, as 
is seen from (flQl) . 



4-2. Omega particles 

Here, we will find the optimal shape of an omega particle, which, according to Section [3l 
will optimally interact with evanescent linearly polarized waves and reactive linearly 
polarized fields. Let us consider a canonical omega particle P HOl E] whose shape is 
shown in Figure O Similarly to the canonical chiral particle considered above, this is 
a connection of two small antennas: an electric dipole and a small loop, but now the 
particle is planar and the induced electric and magnetic moments are orthogonal. 




Figure 2. Canonical omega particle. A uniaxial particle with the polarizabilities 
defined in Table [T] can be formed by adding a second identical particle rotated by 90° 
around the axis z. 



(25) 



For simplification of derivations, we suppose that this particle is excited by a linearly 
polarized evanescent plane wave with the fields 

Et = EyYo. 

In this case, the incident electric field is parallel to the electric dipole antenna and 
the magnetic field is orthogonal to the loop plane. Because the only difference with the 
earlier considered canonical chiral particle is the mutual orientation of the wire and loop 
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antennas, we can use the same expressions for the polarizabihties as in f llTp (neglecting 
the electric polarization of the loop) [8]): 

/2 



-J 



^Ztot 



-JOiccVO- 



(26) 
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The fields exponentially decay along the positive direction of the ^-axis and the 
transverse wave vector is directed along x, with kt > fco- The wave impedance for 
this evanescent wave is Zte = i^/^o//3, where 



(3 = nko = jn"ko = jko, 



1-1 

\ki 



(27) 



is the imaginary propagation constant of the evanescent wave, and the imaginary value 
of n = jn" represents the effective refractive index. From f|T2l) we can write the power 
absorbed by the particle as 



P 



i|Et|2lm|a^: + 2mm 
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Using the expressions for the polarizabihties of omega particles ( 126|1 . we are ready to 
find the absorbed power at the particle resonance: 



P 



2' ^ R 



tot 



(29) 



For the resonant particle the stretched wire length is approximately A/2, meaning that 
the dipole length and the loop radius should satisfy the relation 

A 



/ + vra 
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Substituting in (129|) . we get 



(30) 
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(31) 



Equating the derivative of P with respect to / to zero, we can easily find the optimal 
values of / and a for which P has the extremum value: 

A 



opt 



^opt 



" + 1 - ^n"{n" + 2) 
n" + ^n"{n" + 2) A 



(32) 



TT 
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From these relations, it is clear that for different values of the transverse wave number 
(different decay rates) the optimal dimensions are different. For the interesting special 
case in which the effective refractive index equals unity in the absolute value: n = ±j 
(corresponding to kt = \/2ko), we note that the optimal "scaling factor" koS/l = 1, 
and the necessary condition for the polarizabilities of the optimal omega particle in (l26ll 
reads 

?7oaee = ±Q = — Omm- (33) 
^0 

This is again the same balance relation as for the optimal chiral particle excited by 
a plane wave in free space (refractive index n = 1), telling that the optimal response 
corresponds to all the normalized polarizabilities being equal. 

5. General definition of optimal and balanced particles 

The above two examples of reciprocal bianisotropic particles, optimized for the 
maximum power extracted from the incident fields (Section H]), show that the optimal 
response is related to the condition of balanced reactions to electric and magnetic fields. 
According to flMl) and fl33|) . all the normalized polarizabilities of these optimal particles 
are equal. Some examples of special properties of particles with balanced response to 
electric and magnetic excitations are known from the literature. For the simplest case 
of the absence of magneto-electric coupling, the balance condition ?7oaee = (l/''7o)«mm 
corresponds to the Huygens particle, also called self-dual particle, e.g. [iHl UHl [201 121] . 
These particles have the interesting property of zero backscattering amplitude. Special 
properties of small chiral particles with balances polarizabilities are also known, e.g. 
[H H]. In addition to the property of maximized or minimized strength of particle- 
field interactions, these particles are most suitable for realizing artificial materials with 
negative refractive index using the "chiral route" [221 [231 \2M [251 [261 [27] . It appears 
that the optimal chiral particle exhibit chiral effects with the extreme strength: Optical 
activity in chiral media is caused by the difference in the propagation constants of the 
circularly polarized modes with the opposite handedness, and composites of balanced 
chiral particles can be engineered so that these two propagation constants differ by 
sign |23l [2^ . In particular, it is possible to realize media with the effective refractive 
index n = ±1 for the two modes (in the assumption of small losses). Furthermore, 
the use of optimal chiral particles makes it possible to design single-layer polarization- 
transformation sheets, e.g. [28]. These examples of balanced chiral particles show that 
the optimization for the maximized power, extracted by the particle from the incident 
field, gives rise to a number of extreme electromagnetic properties of such particles and 
their ensembles. 

In this paper we have identified the classes of optimal particles for the most general 
uniaxial bianisotropic particles, including nonreciprocal particles. We have seen that 
the requirement for optimal interaction of omega particles with exciting fields also leads 
to the condition of balanced polarizabilities for the special case of the effective refractive 
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index equal to n = ±j (unity in the absolute value). We expect that not only the known 
chiral particles with balanced polarizabililies show extreme and useful electromagnetic 
properties, but also particles of all the other classes, provided that their polarizabilities 
are balanced. This expectation is supported by our recent studies of scattering by 
individual uniaxial particles and infinite periodical planar arrays [291 130] . In these 
studies, we have found that conditions for zero forward or backward scattering from 
general uniaxial particles lead to the requirement of the balanced polarizabilities. Also, 
it appears that the requirements of perfect absorption, full transmission or full reflection 
from arrays of small particles (both reciprocal and nonreciprocal) again leads to the 
same condition of optimal balanced polarizabilities. Based on the above observations, 
we suggest to introduce the notion of the balanced particle as the particle with balanced 
polarizabilities for all classes of uniaxial particles, listed in Table [H This condition can 
be expressed in the general 

rjOttee = ittem = ± — ttmrn- (34) 

Vo 

The magneto-electric parameter can be any of the four possible coupling parameters: 
the chirality parameter k, the omega parameter Q, the Tellegen parameter x cind the 
velocity parameter V. From the known extreme properties of balanced chiral particles in 
the field of propagating circularly polarized waves, we expect that the balanced particles 
of other classes will show extreme responses when excited by the incident fields of 
other polarizations and by reactive fields, accordingly to the particle classes identified 
in Section [3l Finally, we note that in many cases the balanced particles satisfy the 
conditions of maximal or minimal absorbed power, meaning that they are also the 
optimal particles for interactions with incident fields of various types, according to the 
classes identified above. 

6. Conclusion 

In summary, we have introduced the concept of balanced and optimal particles and 
composite materials for interactions with propagating or evanescent fields of linear 
and circular polarizations. The optimal particles extract maximum power from the 
incident electromagnetic fields (for a given overall size and the resonant frequency of of 
the inclusion). Apparently, these particles are also the optimal radiators of power. 
We have shown that the most interesting optimal reciprocal particles have all the 
polarizabilities equal (or they differ only by sign). We call such sets of polarizabilities 
"balanced" . It is expected that the use of balanced and optimal particles especially in the 
design of artificial electromagnetic materials will allow optimization of electromagnetic 
performance of various devices, where the material or layer is interacting with various 
types of electromagnetic waves (antennas, absorbers, sensors, lenses, and so on). 
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